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Abstract 

We calculate the dissipation of the baryon number after the electroweak phase transition 
due to thermal fluctuations above the sphaleron barrier. We consider not only the classical 
Boltzmann factor but also fermionic and bosonic one-loop contributions. We find that both 
bosonic and especially fermionic fluctuations can considerably suppress the transition rate. 
Assuming the Langer-Affleck formalism for this rate, the condition that an initial baryon 
asymmetry must not be washed out by sphaleron transitions leads, in the Minimal Standard 
Model (sin 9w = 0), to an upper bound for the Higgs mass in the range 60 to 75 GeV. 
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1 Introduction 


The question about the origin of the baryon asymmetry of the Universe (BAU) has 
recently gained much interest. Many different models of how the BAU was created are 
being discussed in the literature (for reviews see e.g. 0), some of them considering 
BAU generation at the GUT stage of the Universe, others favoring the generation 
during the electroweak phase transition. 

Whatever the mechanism was which led to the BAU at early times, the resulting 
asymmetry might have been eliminated by baryon number violating processes in the 
electroweak theory after the phase transition. Such processes are possible due to the 
anomaly of baryon and lepton currents 0 and the non-trivial topological structure 
of the Yang-Mills theory. This feature was discovered in 1976 by Faddeev 0 and 
Jackiw and Rebbi 0], who found that the potential energy is periodic in a certain 
functional of the fields, the Chern-Simons number Acs- Topologically distinct vacua 
of the theory are enumerated by integer Nqs- In the electroweak theory those vacua 
are separated by an energy barrier whose height is of the order of rri\y / a where m\y 
is the mass of the W boson and a = g 2 /(47t) is the 577(2) gauge coupling constant. 

Transitions from one vacuum to a topologically distinct one over this barrier 
change the baryon and lepton number by one unit per fermion generation due to 
the anomaly of the corresponding currents. If we assume in accordance with the 
standard model that B — L (baryon minus lepton number) is conserved and that 
there is no primordial excess of say, antileptons, then these transitions can cause the 
BAU erasure as mentioned above. Hence it is necessary to know the transition rate 
of such processes. While the baryon number of the Universe today (Bo) is about 10 -9 
to 1CU 10 (relative to the number of relict photons), in models generating the BAU at 
the electroweak phase transition the number of produced baryons (Bj J per photon 
is of order 1CT 5 [i|, ||. (The precise value is not important for our calculation, see 
below). Thus the ratio Bq/Bt c which describes the dissipation of the BAU should 
not be significantly lower than about 1CU 5 , otherwise the initial baryon excess is not 
large enough to explain the present day BAU. 

In principle to obtain the value of this ratio one has to integrate the rate of the 
baryon number violating processes over the temperature from T = T c to T = 0. In 
practice, however, the rate is very strongly suppressed at ordinary temperatures 0 
so that only a short range below T c contributes to the erasure of the BAU. While 
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at low T the rate is dominated by tunnelling processes, at higher temperatures the 
energy barrier can be overcome thanks to thermal fluctuations 0 B 0 - This thermal 
transition rate can be evaluated by the semi-classical formalism of Langer |)j] and 
Affleck pi . 

A key role in this calculation is played by the static classical field configuration 
which corresponds to the top of the energy barrier, having Chern-Simons number 
iVcs = This configuration was first found by Dashen, Hasslacher and Neveu JTT 


and rediscovered in the context of electroweak theory by Klinkhamer and Manton, it is 
called sphaleron |l2j. Its energy E r i ass enters the transition rate y(T) = A(T ) e~ Eclass / T 
via the classical Boltzmann factor and is usually the dominant contribution to 7 , 
i.e. in most cases |lnA(T)| < Ed^/T. The prefactor A(T ) contains contributions 
coming from fermion and boson quantum fluctuations about the sphaleron. In || 
the rate was calculated considering the classical and zero-mode contributions, while 
the determinant of non-zero boson fluctuation modes and the fermion determinant 
were set to unity. The result for 7 was so large that any initial baryon excess would 
have been washed out by sphaleron transitions after the phase transition. Therefore 
it is of interest whether quantum loop corrections could help to preserve the baryon 
asymmetry. 

Several investigations which consider loop corrections have already been made. 
Bochkarev and Shaposhnikov [;5j. |j] included the boson fluctuations through an effec¬ 
tive potential of the Higgs field. They obtained that the transition rate is sufficiently 
suppressed if the Higgs mass is below an upper limit of 45 to 55 GeV. A direct com¬ 
putation of the bosonic determinant over non-zero modes was made in JT3| by using 
an approximation technique [|T|], exact calculations were performed in |H| |T 6 j . 

All these calculations were based on the high temperature limit in which the four¬ 
dimensional fluctuation matrix can be replaced by the three-dimensional one and 
fermions decouple completely. Although parametrically this limit is reasonable it 
need not necessarily be justified numerically. In this paper we go beyond the high 
temperature approximation which corresponds to taking into account the fermion 
determinant (suppressed in the formal high temperature limit) and to calculating the 
full four-dimensional bosonic determinant. We generalize the preceeding calculations 
to arbitrary temperatures. First, we include fermion loops which previously have 
been altogether neglected, second, we evaluate the fluctuation determinants for finite 
temperature considering the full sum over Matsubara frequencies. Both contributions 
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seem to be quite essential numerically. 

The fermion determinant has been considered in detail in our previous publication 
B the bosonic one is the aim of this work (see also ||I8|). We find that the result is 
significantly influenced by terms which vanish in the formal high T limit, especially by 
the contribution of the fermion fluctuations. Actually, the polarization of the Dirac 
sea of fermions in the classical sphaleron background held adds up to about 30 % to 
the sphaleron energy. Therefore, the fermion determinant which was put to unity in 
o, m m leads to a strong additional suppression of the transition rate. 

Let us remark that in an abelian (l+l)-dimensional model the transition rate can 
be calculated analytically. This has been done in |0| for the boson and in [0] for the 
fermion loop correction. In (3+1) dimensions, however, one has to resort to numerical 
methods. 

For the evaluation of the boson determinant we use the same method as for 

It is based on the computation of the complete 


the fermion determinant in 17 


(discretized) spectrum of the fluctuation operators. All the relevant quantities can 
subsequently be calculated for any temperature T by suitable summations over the 
eigen-energies. Including the loop corrections into the formula for the transition rate, 
we finally obtain the ratio B 0 /Bt c which is a measure for the erasure of the BAU. 
We find that both bosonic and fermionic fluctuations suppress the rate considerably, 
especially for a low mass of the Higgs boson and a large top quark mass. For a top 
quark mass in the range m t = 150 to 200 GeV, in accordance to recent experimental 
results |2T| , the condition that the BAU must not be washed out by sphaleron transi¬ 
tions leads within the framework of our one-loop calculation and the Langer-Affleck 
formalism to an upper limit for m# in the range between about 60 and 75 GeV. 

Another goal of this work is the recalculation of the boson fluctuation determinant 
in the high T limit since the results of the two existing calculations [0, 0, based on 
different analytical and numerical techniques, deviate from each other. Although they 
show the same tendency for low Higgs masses, no satisfactory quantitative agreement 
was found. Our numerical method is based on the diagonalization of the fluctuation 
operator as described above; it differs significantly from those used in 115, [TTitl so that 
our study can be considered as independent. We find that our results for the boson 
determinant in the high T limit agree with the results of |]TO| up to about 10% while 
there is a larger deviation from the ones of 0- 

The paper is organized as follows: In Section 2 we introduce into the notations 
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and conventions of the model, in Section 3 we apply the Langer-Afheck formalism to 
the baryon number violation processes. The renormalization is discussed in Section 4, 
followed by the treatment of the temperature dependent parts of the fluctuations. The 
evaluation of the baryon erasure B 0 /Br c is done in Section 6. Numerical results and 
checks for the computation of the fluctuation determinants are presented in Section 
7. In Section 8 we proceed to the evaluation of the sphaleron transition rate and 
deduce the upper bound for the Higgs mass. We also investigate the applicability 
of the framework of our calculation. Finally we summarize the results and draw our 
conclusions in Section 9. Technical details about the computation of the discretized 
spectrum and the spectral densities are treated in the appendices. 


2 The model and parameters 


We consider the minimal version of the standard electroweak theory with one Higgs 
doublet which is Yukawa coupled to left handed fermion doublets and to right handed 
singlets; in the following we write only one doublet and one pair of singlets for brevity. 

We shall work in the limit of the vanishing Weinberg angle, i.e. the theory is 
reduced to the pure SU(2 ) case (without the 77(1)). This idealization does not seem 
to be significant [^]. The Lagrangian is thus 

c = - 

+ ( 2 . 1 ) 


with the covariant derivative D M — iA M , A ^ = |A“r a , and the held strength 
= \F^r a = i[D^D v ] , F^ = d,Al - d„A* + e^A^Al. M is a 2 x 2 matrix 
built of the Higgs held components $ = (^o) and the Yukawa couplings h u , h^; 


M 


h u $°* h d § + \ 
-h u $+* h d $° ) ' 


( 2 . 2 ) 


VtL means the SU(2 ) fermion doublet 

^L = 5 ( 1 - 75 )^ = 

and with i/j r we denote the pair of the singlets 

4 >r = 5(1 +75)V> = 



(2.3) 


(2.4) 
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The masses generated by the non-vanishing vacuum expectation value (0|$|0) = 


:G) are 


gv h u4 v 

m w = — , m U)d = -~j=- , and m H = Xv 


for the gauge boson, the fermions, and the Higgs boson. 

We prefer to work in terms of dimesionless rescaled quantities 

mw 


—> m^x 11 


A a —>• m w A a , 


<f> 


$ ; 


v /2 g ’ 


using the following representation of the Dirac matrices 


(2.5) 


( 2 . 6 ) 


7° = 


0 1 
1 0 


7 — 


0 Oi 
-Oi 0 


7 5 = 


-1 0 

0 1 


the fermion Dirac spinors can be reduced to two components: 


V’l -> m% 2 j , 


if)R —> m 


3/2 | 
W | 


f 0 N 

<$>R, 


In this representation the Lagrangian (|2.1|) is 


C = 


+ - (TiD^tjjL + libido + (Jidi)i; R - i^ L Mii R - 


with the mass matrix 


M = 


2m 


w 


m u $°* mrfd> + 
—m„$ + * md&° 


and 


v H = 


m w 


(2.7) 


( 2 . 8 ) 


+ i(£>„$)*(£>"$) - ii£(*** - 4) 2 ) (2.9) 


( 2 . 10 ) 


3 The baryon number violation rate 

As it is well known || the baryon and lepton numbers B and L are not conserved in the 
standard model of electroweak interactions due to the anomaly of the corresponding 
currents jg and jjh 

N 

d,f B = d„j£ = ( 3 . 1 ) 

Here N g = 3 is the number of fermion generations. We dropped the contribution of 
the 17(1) gauge field to the anomaly since we work in the approximation of vanishing 
Weinberg angle sin 6w = 0. 
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Integrating this anomaly equation one finds that in any process the change of the 
baryon and lepton numbers is related to the change of the Chern-Simons number by 


AB = AL = N g AN CS , 


(3.2) 


where the Chern-Simons number Nqs is defined by 

Ncs = f d 3 r Jl 0 (r), 


(3.3) 


iW = 


167T 2 




d K 1 ' = 

11 64vr 2 


(a&a; + , 

1 


pvpa pa pa 
1 p J v ± pa 


(3.4) 

(3.5) 


The semiclassical description of the processes with fermion and lepton number viola¬ 
tion is based on the existence of an infinite number of classical vacuum configurations 
labeled by integer values of Nqs- These vacua are separated by potential barriers 
which can be overcome either by quantum tunnelling or by real processes in the 
Minkowskian time due to thermal fluctuations. We shall be interested in tempera¬ 
tures at which the real-time transitions dominate. The rate of the thermal transitions 
at temperature T between adjacent vacua is roughly given by the Boltzmann factor 
exp(—A c i ass /T), where E c i ass is the energy of the sphaleron, which is the field config¬ 
uration at the top of the barrier between the two vacua. The transition rate T with 
the preexponential factor is given by the Langer-Affleck formula 
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r = 


|cn_| Im Z. 


sphal 


7T 


(3.6) 


Here Z 0 is the partition function computed in the semiclassical approximation around 
the vacuum, and Z spha i is the partition function obtained by semiclassical expansion 
about the sphaleron solution. Since the sphaleron solution is a saddle-like point of the 
potential energy functional, the quadratic form of fluctuations about the sphaleron 
has a negative mode uA < 0 so that the formal semiclassical expansion about the non- 
stable static solution gives a complex contribution to the partition function Z sphal . 

In the Weinberg-Salam model the sphaleron solution in the temporal gauge A 0 = 0 
can be found as a stationary point of the energy functional 


A'class 


rn w 


(TV 


i(f)“) 2 + 1(A*)'(A*) + i^(4> f 4> - 4) 


32 H v 


(3.7) 
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which leads to the classical equations 


(D,F tl ) a - |($V(A$) - (A$)V$) = 0 

(A 2 - X(^-4))$ = 0 . 

The sphaleron solution is assumed to have the following form (hedgehog) 

1 


(3.8) 


A i( r ) = e aijnj 


A(r) , B{r) C{r ) 

\^ai W'aV'i) 


<f>(r) = 2 H{r ) + iG(r ) n ■ 0 


(3.9) 


ffere r = rn, 0 are Pauli matrices, and the profile functions j\,B ,C ,G ,11 can be found 
numerically by solving the classical equations. 

The spherical symmetry of this static solution is preserved under time-independent 
gauge transformations of the form 


f/(r) = exp iP{r ) n • 0 


(3.10) 


with an arbitrary function P(r). One of the five profile functions could be completely 
eliminated by using this gauge freedom, but since in this case the remaining functions 
are not necessarily regular at the origin and at infinity, which is required by our 
numerics, we use all five functions. 

It should be mentioned that the expression for the Chern-Simons Number given in 


eq. (|3.3|) is not gauge invariant, so it is only well-defined if we require the fields to be 
continuous at infinity. In this case Nqs is determined up to an integer number which 
is the winding number of a possible gauge transformation. The final results should be 
independent of the choice of the gauge. Since we did not exploit the gauge freedom 
to eliminate one of the five profile functions, we can verify this gauge invariance 
numerically. This provides a powerful non-trivial check of our performance. 

We can rewrite the general Langer-Affleck formula eq. (p.6|) in the form 


T = 


7Tl\yLU _ 

2 t r 


j VaVipV^Vip exp (~S[A + ga,$ + ) 


exp 


(-S[A (0) + ga, <f> (0) + g(p, -0] 


2 nd 

order 


(3.11) 


ffere the Euclidean action 

1 


I r/3m w r r 

S[A,4,y,V<] = - / dt d 3 r i(F^) 2 + ±(.D„4) , (rV4>) 

(J J 0 J 


g f 


+ M I 4($ t $-4) 2 + A (0l Xr) 


D 0 + idiDi 
Mt d 0 


M 

- mdi 



(3.12) 














is expanded to second order in fluctuations around the sphaleron configuration A, 
in the numerator and around the vacuum configuration db°£ in the denominator. 
(3 — 1/T is the inverse of the temperature T. 

In zeroth order (no fluctuations) S reduces to (3 times the classical energy ( |3.7|) of 
the sphaleron or the vacuum configuration, hence the transition rate V contains the 


Boltzmann factor exp[— (3(E c i ass — mentioned above, where we have E^ s = 

E dass [A(°\$M}= 0 . 

In fact, formula (|3.11|) should be modified to take into account the gauge fixing, 
renormalization and special treatment of zero modes. Let us start with the gauge 
fixing. Following |H| we shall work in the background R g gauge defined by adding 
the term 

5 ((4<v)“ + J(*VV - A“*)) 2 (3.13) 

to the quadratic form for the fluctuations a, <p. D is the covariant derivative with the 
background field A. In this gauge the Faddeev-Popov determinant is 


do) 


k FP — 


det(—+ /Cfp) | with 


«fp = - A 2 + j. 


(3.14) 


The quadratic form of the action expanded about the sphaleron solution in this gauge 

is 


S^S — ~J d A x a“(—- D\ + ;| < f )t< f ) ) a 




i<((A$) t rV - Ar“(A$)) + yf ($V + 


- —(&r a ip - Ar a <F) 2 + <pU-d% - D] + _ 4 A y (3.i5) 

16 ' S / 


2 (A A 


f j 'r) 


d 0 + i(7 j Dj 

Aft 



where M is the mass matrix ( 2.10|) with the background Higgs field $. The integration 
over a 0 yields the inverse square-root of the Faddeev-Popov determinant, 1/afp- 
The spherical symmetry of the sphaleron solution leads to certain symmetries of the 
quadratic form of the action. In order to make these symmetries explicit in the 
Higgs sector, we prefer to work with the complex Higgs doublets in terms of four real 
components. Any complex doublet £ can be represented in the form 

^£2 + *£i 


£ = 


.£4 — *£■? 


= £ r + 
W 


'(P 

. 1 / 


with 


r M = ( ± * 0 ’ 1 ) (h- = P---, 4 ) ; (3-16) 
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where is a real four vector. One has 


e C + = 2LC 


M A* 


i (fVC - CV“£) = 2r ] a ^ fl ( u 


where r}^ v is the’t Hooft symbol 


^=4i‘ r 


TIT T 
\ fj, 1 / 





(3.17) 

vtj = «»« • 

(3.18) 

- (D>.i)p* (; 

J can be 

1 a A a 

(3.19) 


written as 

(-Da^) m = {Dxj^v with {Dx)fj,u = ^,wd) 

Joining the fluctuations of the gauge and Higgs helds into a 13 component vector 
(a“ 99 ^) we can rewrite the quadratic form of the action (|3.16|) in the form 


5 {2) S= [ 
J 0 


/ 3m w 


where 


d 4 x 


b(a <f)(~do + ^bos)l 


dfij 


+ {^L ^R)(d 0 + JJferm) j , 


(3.20) 


^hos 


C a f> W a 

27 IJ r Y IV 

Wl Hr 


_abc t^c 
r ij ? 


£“/ = Jp(-(^) a& + 1 <5 a6 i> a l> a ) + 2e c 

'Hnv = -(Dl)^ - i(l - l4+ i V((! + ~ 

Wl = C(Ai>)a , 


J^ferm — 


mDi \v F ^^r+ 




Tcq<9; 


(3.21) 


In order to preserve the spherical symmetry of the fermionic Hamiltonian 7df erm we 
have taken equal masses for up and down fermions: 

m u m d 


m u = m d , 


u F — 


(3.22) 


7Jl\v m\y 

The physical significance of this approximation will be discussed later. 

Using the quadratic form of the action Q3.20| ) we can perform the functional inte¬ 
gral in the Langer-Affleck formula for the transition rate (|3.11|) and get f| H 

,2 T / _,J0) Kr 

(3.23) 


__ 87 T V 3 3 AAfP ^bos ^ferrri r_/opbare] 

^ 9^ (g 2 / 3) 3 ^ tr ^ rot ( ' n '' - eX Pl- T-^classJ • 


4°p Khos 41 
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Here ^ is the classical energy ( |3.7p of the sphaleron solution with bare (unrenor- 
malized) parameters. The (dimensionful) quantity V is the physical space volume 
arising from the integration over the three translational zero modes of the sphaleron. 
Moreover the sphaleron has three rotational zero modes. The effect of all zero modes 
is taken into account in eq. (|3.23|) through the factor [|l3j iV t 3 r A^ r 3 ot with the 

Jacobians 


N tr = 


N rn t — 


67T 

6n 


d 3 r («) 2 + 

d 3 r ( (rHj, - r }ri )(F?jFS + (D,9) 3 (D,i)) - e m A K 


(3.24) 




where A is the solution of the equation 

(-A 2 + . 


(3.25) 


The determinants kfp, Kbos, «ferm in eq. (|3.23|) arise from the Gaussian integration 
over fluctuations with periodic time boundary conditions for bosons and antiperiodic 
boundary conditions for fermions. We can write them in the following form 


Kpp — 


0 _,fpA 


det(—<9 q + /C F p) 2 = n( 2 sinh(|m^w" )J 

n 

det'(—9g + /C bos ) 2 = R (2 sinh ( ^m w uj n 

n 

^ferm (1 ('t (()() T Tfferm) COSh( - 


^bos 


(3.26) 

(3.27) 

(3.28) 


where od,u; p , and e n are the eigenvalues of /Cbos, /Cpp, and 7df erm , respectively. Our 
numerics is based on a calculation of the eigenvalues of discretized versions of these 
operators; details can be found in Appendix A and in |T7|. The primed product in 
eq. (|3.27|) means that the zero modes are omitted, the negative mode however, 
contributes to ftbos- 


4 Renormalization 

Equation (|3.23|) for the transition rate contains ultraviolet divergences arising from the 


infinite products in eqs. ( p.26| - [3728|) . These divergences are removed by the renormal¬ 
ization, where it is sufficient to renormalize the theory at zero temperature. Keeping 
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this in mind we split the right hand side of eq. (|3.23|) into the temperature dependent 
finite part and the divergent part corresponding to zero temperature. We write 

8tt 2 V 


r = 


fnw\uj -1 


47 rsin(fm w |u;_|) g 6 f3 3 


(iVniVrot) 3 exp [~p(E, 


bare 

class 


E\ 


-T=0 

FP 


+ £‘T(T) + ££° + Eg>(T) + Ef~° + E£™(T)) 


where the ultraviolet divergent T — 0 terms are 


7lT=0 

7 FP 


■\m w 


(5>r-E 


CO, 


ppO 


n n 

E141) • 


f = -W(E 




The temperature dependent terms 
B‘*»p ( r) = -i(5>(i 

n 

crm = +j(e" mi- 

n 

Kz r m = -i(E >“(!+■ 


e ln(l — e 


^—f 3 m w un^ 


,-fSm w |e„p 


-EMi - 

n 

E ln (! + 


-0m w u° 


D -0m w |e°l 


(4.1) 


(4.2) 


(4.3) 


are finite and vanish at zero temperature. Here XT stands for the sum over all non¬ 
zero, non-negative modes. 

We shall perform the renormalization of the zero-temperature contribution using 
the proper-time representation for the quantities (|4.2|) with the cutoff parameter A: 

e fp°( A ) = J A _ 2 ^ Tr ( exp [ t)C,pp} - exp[-t)C <E }) , (4.4) 

E lTm( A ) = J a _ 2 ^2 Tr ( ex P[- m ferJ - expH^iE) 2 ]) • (4.5) 

In the case of E£=° the proper time representation is modified to suppress the negative 
mode contribution: 

£bof(A) = J a _ 2 ^jr 2 { Tr (exp [-t/C bos ] - exp[-t/c£i]) + (1 - e _faJ -)} (4.6) 
In the limit A —> oo these integrals diverge since for t —> 0 

Tr (exp {-tK] - exp[-f/C (0) ]) = at~ x ' 2 + bt 1 ' 2 + • • • , (4.7) 
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where /C can stand for /Cbos, /Cfp, /Cf erm = Aff erm . We write the divergent pieces of the 
r.h.s. of eqs. i|l. 1| |l.(i|i as 


pdiv _ fffw_ [ v ™ dt 

FP 40F Ja-2 A 2 


— [Tr exp(—f/C FP )] div , 


77 'div _ 

-^bos 


77 'div _ 

-^ferm 


mw 


dt 


Ay/ii J a - 2 f 3 / 2 


Tr exp(—f/C bos )] div , 


mw r™* dt / +jr xi 

y A _ 2 £ 3/2 t Tr exp ( A/Cferm)]div , 


where we define 


[Tr exp(—t/C)] div = at A 2 + , 


(4.8) 


(4.9) 


and v Teo mw is the renormalization scale which is determined below from the value of 


the Higgs pole mass. Performing the small t expansion (|4.7| ) and integrating over t 
in O we find 
m w 


77 'div 

U/pp 


^bt(A) = 


(A) = 


64(/r 2 


d 6 r 


-fz/ 2 

2 V ren 


A 2 )($t$ _4) 


hi (£) B<£) 2 +)h 4, ' I> - 4 > 2 +K*’* - 4 >) 


mw 

64(/r 2 


d 3 r 


-(V 2 

2 \ ren 


A 2 )(4 + >/£)(*'* - 4) + In (£) (f (f“) 


+ 6(A't) , (Bi*) + 4(4 + vl + 4)(* f * - 4) 2 + f(8 + + i4)(it# - 4) 


16 

j3 


Bf “ (A) = S ? / <ft '[ 4!/2 ( !/ L - A 2 )(#t« - 4) + In (£) (i(F“) 2 

+ 2zz 2 (A$) t (A$) + |^($ t $-4) 2 + 4^($ t $-4)) . (4.10) 

The fermionic part in the above formulas is written for one fermion doublet. As 
it was mentioned above, to preserve the spherical symmetry of the Dirac equation 
one has to consider the case of equal masses for up and down fermions. This approx¬ 
imation is justified for doublets where both fermions are light. However, in the case 
of the (t, b ) doublet the approximation m & -C mw -C wi( is more reasonable. In this 
limit the correction from the top quark is half that of the doublet with both masses 


equal to m t |17[- Therefore, in our numerical estimates we take 9 + | massless fermion 
doublets and | massive doublets with mass mt, taking into account that the quark 
contribution is enhanced by three colours. 

We see that the quadratically (A 2 ) and logarithmically (In A 2 ) divergent terms 
are exactly those entering the classical energy functional (|3T7|) . Therefore, they 
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can be combined with the bare constants of the correspondent terms in the classical 
energy — to produce the renormalized constants at the scale of u ren mw • We call the 
classical energy with renormalized constants -Edass- 

What is left after the renomalization is ultraviolet finite, and one can safely put 
the ultraviolet cutoff A to infinity. We call these pieces renormalized energies, 


with 


77>conv 

-^FP 


(A) = 


rn w 

4 


£fp = lim ££°p nv (A), 

A—>00 


f°° dt 


(4.11) 


/ A _2 (3/2 Tr (^32p[ l^Fp] - exp[—//Cpp]) 

r ^ren dt 

-jyj [ Tr exp(-«Kpp)]<uv 


(4.12) 


and similarly for E^ m and E^. I n the case of E^l one should take into account the 
subtraction of the negative mode contribution in (|4.6|) so that 


^bT(A) = 


m w 
4 y/n 


r°° dt 
/a-2 W 2 L 


Tr ( 


exp[ 


-tic 


bos I 


exp[ 


-tic 


(o) 

bos 


+ (i 


— tuj z 


ren dt 

L-* W 2 


Tr exp(-f/C 


bos/ldiv 


(4.13) 


Performing the described renormalization procedure we arrive at the following 
expression for the transition rate ([4.1|) 


7 = v = 


27Tm w \u-\(N tr N rot ) s 
g 6 /3 3 sin(|mw|u;_|) 


exp 


A(£class + + E%? P (T) 


+E, 


ren 

bos 


+ <™ p (r) + ■^ferm + KrZ(T)) 


(4.14) 


Finally we fix the renormalization point u ren in eqs. (| 4 . 8 | , |4 . 1 2 | , | 4 . 13 |) so that the 
renormalized parameter n# coincides with the physical pole mass. In order to ob¬ 
tain the pole mass we have to evaluate the propagator of the Higgs particle in one- 
loop order. Its classical part (in Euclidean space) is given by G~ l {jp 2 ) = p 2 + u 2 H 
(in units of m 2 v ), the fermionic and bosonic one-loop corrections can be written as 
a N c v? Cj) and a u 2 en F hos (u 2 en , p 2 , u 2 H ) with some functions F ferm and F bos 

which are finite in the limit of infinite u t and z/ ren , respectively. 

In most parts of our numerical computation we work with a top quark mass 
substantially larger than mu- From the classical part of G^ijp 2 ) we know that p 2 is 
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of the order of —z/#, while we find (see below) zy 2 en ~ z/ 2 . Therefore, the dominating 
parts of the one loop contribution are terms of the order a N c uf and a N c v 2 stemming 
from the fermionic correction, so for m t significantly larger than m h it is a good 
approximation to drop the bosonic contribution completely (its leading term is only 
of the order a z/ 2 en ~ olv 2 ) and to expand Ff erm in p 2 /z / 2 up to the order 0(p 2 /v 2 ). 
Within this approximation we can also neglect corrections to the pole mass of the 
W-boson. Exceptions from this procedure will be treated below. 

In order to get the Higgs propagator we put the gauge held to zero and obtain 
the relevant (four-dimensional) action: 

S — So + Si 00 p , (4-15) 


with 


»SV) = 


r 


drx 


±(<9^)^$) + §($ f $) 2 - 


(4.16) 


1 


^loop 2 Tr lo >gD*D + S t 


where we used 


-l r - Tr e~ tD]D 

counter ~ ° 


dt 

2 A- 2 t 


rU ren dt ( CL 


A^/tt J a - 2 t \ t 


- t + 6 , 


D ] D = -d 2 + iU^odo + 7 ^)(M±±p + M ti ^)l + M f M, (4.17) 


and the counterterms 


a = 




6= ^ 


647r 3 / 2 


(4.18) 

The vacuum expectation value of the Higgs held in one-loop order is obtained by 
setting $1$ = v 2 = const in the action ( |4.15|) and minimizing it with respect to v. 
The result up to order g 2 is 


g 2 N c 


1,2 = y2 ° + (AAn - ^ 4 (l - C - In £-)) , 


(4.19) 


H 


where v'q = 4 is the vev. on the classical level and C ~ 0.577 the Euler constant. 
Choosing the unitary gauge, we can substitute 


4> = 


' 0 N 
K v + gr}j 


(4.20) 
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The propagator is then given by 


S 2 S 


Srj(xi) 5r)(x 2 ) 

and we obtain up to order g 2 


/ tSi 

rj=0 J yZ7T) 


(4.21) 


n-l( 24 2 I 2 I 9 2 r 

G (P ) = P + U H + 8?r2 

with the function 


n 2 

ren 


l-^(C + In#)^(2 + |4)/(-4) 


(4.22) 


/(x) = b f d/3 In 11 — /3(1 — /3)x| = —1 + Re ( y-artanh J —-—- I . (4.23) 


x 


x — 4 


As mentioned before, we expand the propagator up to terms of 0(p 2 /v 2 ). We obtain 


G-\p 2 )=p 2 + v 2 H + 


g 1 N c u t 
8 tt 2 


4 r 


_ 1 _ 1 p 
1 4 ^ 


(I + c + ln i) 


(4.24) 


Our aim is to fix z/ ren such that the pole mass u p , defined by G 1 (p 2 = —v 2 ) = 0, 
coincides with uh- Hence we have to solve the equation 


'4n = 4 


^(| + c+ ,„ 




uz 


(4.25) 


which determines the renormalization constant z/ ren for given values of z/ t and 
u t > v H - As anticipated, we find z/^ en ~ u 2 . 

Analogous equations follow from the full (non-expanded) propagator ( |4.22| ) and 
from the propagator which includes both fermionic and bosonic fluctuations. We have 
checked that for u t > v H the results for z/ ren obtained with those equations are very 
close to the solutions of eq. ( |4.25| ), so that in this case the restriction to the dominating 
parts up to Q(aN c u 2 ) is a very good approximation, and for our choice of u ren the 
renormalized parameter vh corresponds to the Higgs pole mass very accurately. 

The situation is slightly worse for the case of very large vh (e.g. mu = 350 


GeV). Here no solution of eq. (|4.25|) can be found, so we choose u ren such that the 
difference \u p — uh\ takes its minimum. The deviation, however, is found to be below 
10%. Moreover, we consider such high Higgs masses only for comparison to our main 
results, so that this problem is actually irrelevant. 
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5 Thermal renormalization and spectral densities 


We next consider the temperature dependent terms in the exponent of ( |4.1|) given by 
(|4.3| ). In order to compute these quantities we introduce spectral densities Q"\E) for 
the continuous parts of the spectra of the operators /C..., such that for any function 
f(x ) that vanishes fast enough for x —> oo we have 

Tr (/(K) - /(£<">)) = V /(w 2 ) + r dE e(E)f(E 2 ) . (5.1) 

a ; 4 - „ “ 0 


dicrete 

levels 


/C bos has = 7 discrete levels (six zero-modes and one negative mode), /Cf er m has 
n ferm _ i c ii scre te zero mode 0 , /Cfp has n FP = 0 discrete modes, which by definition 
are not included in the spectral densities Q 

We rewrite eq. (|4.3|) with the help of the spectral densities: 


/*oo 

= ±± jf dEQ hos,FP (E) ln(l — e ~ l3mwE ) , 

/•OO 

E^{T) = -\\ dE Q ierm (E) ln(l + e - pmwE ) - j n 

J 0 


ferm 

D 


In 2 


(5.2) 


where the signs are chosen according to those of E temp (T) given in eq. ( [4.3| ). 

At high temperatures eq. ( |5. 2|) is dominated by spectral densities at large energies. 
It is shown in Appendix B that asymptotically all three spectral densities £? FP , Q hos 
and £? ferm approach constant values (after subtraction of the spectral densities of free 


operators, which is implied in the definition ( |5.1| )), 

lim Q(E) = . 


E^o o 


Therefore, at high T the quantities E temp (T) have a T 2 behaviour: 

/•oo 7J-2 q - 1 1 '2 

£large( T ) = ± 1 / dE g ln M ± e ~Pm W E ) = ± ^oo /_]_ ± 3) 

p Jo 24 mw 

whereas the remaining part of E temp (T) can be written as (see eq. (|B.13|)) 


(5.3) 


=4 


dE ( Q hos ’ FF (E ) - ^ os ’ FP ) ln(l 


0 —/3m w E\ 


= ±- 


/•OO . v 

dE (> os ’ FP (E) - ^ s ’ FP ) In 


]_ _ e ~pm w E 


f3mw 


T^n h ° s ' FP HPm w ) 


££$ P =F Jr#-* 1 * ln(/W) 

/•oo 


T/i 


'small 

ferm 


/•oo x 

-i / dE (V erm (E) - 4 rm ) ln(l + e 

•Z 0 


—(3m w E 


) — —Tl 
) 8 li 


In 1 ™ In 2, 


(5.4) 

(5.5) 
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With these definitions the E sma]l are finite in the high temperature limit (see App. B). 

We notice that at T — 0(mw/g) the local functionals A large (T) can be of the 
same order as the (renormalized) classical zero-temperature energy of the sphaleron, 
therefore in that range of temperature one has to find a new sphaleron solution which 
is a saddle point of the temperature dependent functional 

E'z,m=E d „,+ 4 t(t)+ Egrix) +ism. (5.6) 


Fortunately this functional has the same form as the original U c iass, but with temper¬ 
ature dependent parameters. Therefore its saddle point will be just a rescaled version 
of the original sphaleron configuration. 

Using the expressions for Q computed in appendix B eqs. ( |B,8| . |B.4| ) we arrive at 
mw 


EZJA^-T] = 


g J 

mw 1 

g 2 J 

mwq 

g 2 


d 6 r 


i(^) 2 + |(A$) f (A$) + - 4) 


4 V 1 2 

2t'2 


g t 


32 m 2 w K3 


{\N c v 2 + < + 3 )($t$_ 4 ) 


d 3 r 


d 3 r 


AA + 4(A$) f (A$) + - V) 


^^(4-4g 2 ) 2 


i(A“) 2 + i(A*) f (A*) + - 4) 


A^(4g- 2 -4)^ 


(5.7) 


where in the last expression we change the integration variable, r —> q l r, and use 
the notations 

A(r) = q~ 1 A(q~ 1 r) , $(r) = q~ 1 $(q~ 1 r) (5.8) 

with (compare e.g. P3|) 


q = q(T) 



T 

-J- c 


2 A v H m w 

g 



(5.9) 


The new temperature dependent sphaleron configuration A q , <f> <? which is the sad¬ 
dle point of Af a ss (T) can be expressed in terms of the old zero temperature solution 

A q ( r) = qA(qr) , i> 9 (r) = q<b(qr) . (5.10) 
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Subtracting the vacuum contribution we find 


ElZ, [/V i"; T] - SS,[^ (0) , * 0(0) ; T\ = 9£ci„,[A «>] ■ (s.n) 


If we now replace A, $ by A q , in all parts of our expression for 7 , we find that we 
get the old results again but with m,w replaced by qmw wherever it appears (except in 
the definition ( |5.9| ) of T c ). The vacuum expectation value of the Higgs field becomes 
temperature dependent; it is given by = 2q(^j and vanishes for T —> T c . 

Tims the final result for the transition rate per volume is 


7 = 


2'Kqm w \uj_\(N tv N mt y 


{qm w (3) n ™ exp 


g 6 /3 3 sm(£qmw\u-\) 

- P(E^‘"(T) + K™"(T) + EI^(T)) 

— T exp —^(qEtfass + Ep erm + Epos 


0<l{Ea„, + E% + E{Z + E’r ^ 


ferm J 


m w ^qm w 


(5.12) 


with the prefactor 


T = 


2ir(qm w ) 8 [3 4 \u;- \ (N tr N rot y 


(yf 6 sin(|gmiy|u;_|) 
and the total fermionic and bosonic one-loop contributions are 


(5.13) 


-^FERM 

-^BOS 


~ rnren . zT’small (rp\ 

9-^ferm ~T~ ^f er m K 1 ) 


qmw 


qEg + Eg^T) 


qmw 


+ qE£ p n + £^ a11 


(T) 


qmw 


(5.14) 


6 Dissipation of the baryon asymmetry 

In this section we express the erasure of the BAU after the electroweak phase tran¬ 
sition through the transition rate y(T) following the considerations of 0 (see also 

0 ). 

Transitions over the sphaleron barrier change the baryon and lepton number by 
A B = A L = N g A Acs, where N g = 3 is the number of generations. 

If there were neither baryons nor leptons initially, the transitions increasing Nqs 
would be compensated by processes decreasing N cs so that creation and annihilation 
of baryons would cancel each other. The situation becomes different, however, if we 
assume the initial existence of baryons and leptons. These particles could have been 
created long time before the electroweak phase transtion in the age of GUT or maybe 


19 











during the phase transition. In this case the chemical potentials /i of the baryons and 
leptons are non-zero which leads to an additional term fJ,Ncs in the classical energy 
functional. In accordance with the Le Chatelier principle the transitions will favor 
the wash-out of any particle or antiparticle excess. 

Our analysis will be restricted to the usual electroweak theory where B — L is 
conserved. Therefore, the study of the erasure of the BAU in the context of sphaleron 
transitions makes sense only under the assumption of the initial condition B = L. 

A sphaleron transition with ANqs = 1 creates one particle per fermion doublet. 
We introduce chemical potentials /if (i = 1,...,9) for the quark doublets and /if 
(i = 1,...,3) for the lepton doublets. Then each transition increases the energy 
by AiVcs (j2i =i h? + XEi l 1 i') > he. the classical energy functional fl3.7|) has to be 
replaced by 


25 


-F'class * -Edass T N Qg Etf + E/* 


9 3 

Q , ,.L 


( 6 . 1 ) 


\i =1 


i= 1 


This fi contribution leads to the asymmetry in the Langer-Affleck formula (|5.12|) with 
respect to transitions increasing and decreasing Nqs- We have to set Ncs = +\ for 
transitions which increase the fermion number and Nqs = — f for transitions which 
decrease it. Since the baryon and lepton densities considered are very small we can 
restrict the transition rate to terms linear in /i 


7 + = 7 


7 - = 7 


i-f (x>?+I>. L 

Z \j=l 2=1 


P 


1 + a (E v? + E vf 

Z \i=\ 2=1 J 


( 6 . 2 ) 


so that the baryon and lepton number dissipation, given by the difference of these 
rates, reads: 

7R rU, /JL „ 1 .\ 

(6.3) 




Now we have to express the chemical potentials through the particle numbers in 
order to get a differential equation for the baryon number decrease. For small /i 
standard Fermi-Dirac statistics yields the relation 


N (V>) = 


(6.4) 
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where N is the number of fermions, / the number of degrees of freedom, m the mass 
of the fermions and 1(a) is given by 


12 r°° dx x 2 — \ a 2 
7 t 2 Ja 1 + e x \Jx 2 — a 2 ’ 


(6.5) 


with the properties 7(0) = 1, 7(oo) = 0. A lepton doublet has three degrees of 
freedom (two for the lepton and one for the neutrino) while there are four for a quark 
doublet. The masses of the leptons and the light quarks are much smaller than the 
critical temperature T c = 1 //3 C , i.e. f3 c m < 1 so that I(j3 c m ) ~ 1. We apply this 
approximation also in the case of the top quark since it has almost no influence on 
the result for the baryon number wash-out. Hence we obtain for the lepton and the 
light quark doublets: 

l _ 2/? 2 T Q _ 3f3 2 „ 

y ^ ’ l^i — 2 y ’ (6.6) 


where L;, Qi are the numbers of leptons and quarks of a fixed doublet i. Substituting 
eq. (|6.6|) into eq. (|6.3|) and using L = B — Q /3 we obtain: 


dB 

dt 


= -7 (T) VN g (3 


3(32 Q + %- l) = 7 (T) N g f3 3 B 


2V 


V 


(6.7) 


Standard cosmology gives a relation between time and temperature 


26, 27 


t = 


45 


m P T~ 2 = CT 


16t r 3 lV(T) 


-2 


( 6 . 8 ) 


where N (T) is some number related to the number of degrees of freedom of the 
therm a h zed particles at the temperature T (for our range of temperature it is usually 
^|i), mp = 1.5 • 10 17 my/ is the Planck mass and hence the constant C is given by 
C ~ 5 • 10 15 m^. Substitution yields 


1 dB 
~B ~dT 


13 NgC 


7CO 

J'G 1 


(6.9) 


which can be integrated to 


B(T) 


B(T C ) exp 
B(T C ) exp 


- 13 NgC MT) qi ( q ) ) 

Tl Jo (1 - q 2 y/ 2 q \ ’ 


( 6 . 10 ) 
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with q(T) = yl — T 2 /T 2 . This is our final result; it describes how the erasure of 
the BAU, measured by Bq/Bt c , can be obtained by an integration of q(T) over the 
temperature. In the next two sections we present numerical results of this ratio from 
which we deduce an upper bound on the Higgs mass. 


7 Numerical Results 


In this section we present the results of our numerical calculations. We take mw = 83 
GeV, g = 0.67 which is the physical value of the coupling constant and vary the top 
quark mass in the range 150 to 200 GeV, i.e. around its recently stated value of 174 
GeV |H]. The only unknown physical parameter left is the Higgs mass rn h■ We 
discussed the fermion loop contribution already in [|17] so in this work we will focus 
our attention on the evaluation of the bosonic loops. 

As we explained, all our results are obtained by a diagonalization of the boson 
fluctuation matrix K\ ms and the Faddeev-Popov matrix /Cpp ( |3.21| ) in a discretized 
basis. The spectrum of the matrices must not depend on the choice of the gauge for 
the classical sphaleron fields. Since we evaluated the matrix for an arbitrary gauge 
with non-vanishing C-field (see eq. (|A.25|) ), it is possible to check the invariance of 
the spectrum under gauge transformations eq. ( |3.1(J|) . We find that the eigenstates 
are indeed invariant under gauge transformations if their energy is less than about 
0.8P max where P max is the numerical parameter which restricts the momentum of 
the eigenstates and renders the basis finite to allow a numerical diagonalization (see 
eq. (|A.11| )). Eigenstates with energies close to P max can be gauge dependent which 
is due to the finite numerical box and should not be encountered in our calculations. 
Hence we always have to choose P max large enough so that all eigenstates which enter 
the calculations have energies less than 0.8 P max and, of course, no result changes if 
-Pmax is further increased. 

Another check of the spectrum consists in an investigation of the negative and zero 
modes. The negative mode appears in the grand-spin K = 0 sector of the fluctuation 
matrix /Cbos- We checked that its energy is gauge invariant and independent of the 
box parameters P max and R. Moreover our results agree with the ones obtained in 
[p8[ , [2fj|| where the negative mode has previously been calculated. The zero modes 
can be identified in the K = 1 sector. Due to the spherical symmetry in this sector 
each state is (2K+l=3)-fold degenerate so that we find two threefold degenerate 
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states with zero eigenvalue. Numerically the modulus of the eigenvalues was found to 
be below 10 -3 which shows that the diagonalization reproduces the zero modes with 
excellent accuracy. The eigenfunctions of the zero modes can be evaluated analytically 
in terms of the sphaleron background fields [|13|]. We compared these functions with 
those which we obtained as zero mode eigenfunctions in the diagonalization and again 
found a very good agreement. 

Beside these investigations of specific eigenstates we checked a property of the 
spectrum as a whole. For low values of the proper time parameter t we consider the 
expansion 

_ e -t{u° n ) 2 ^ = Tr (exp[—f/Cbos] - expf-i/Cj®]) 

= ab os t~ 1/2 + b hos t 1/2 + c hos t 3/2 H- (7.1) 


EG 


-tut 


where the coefficients are given in eqs. ( [B. 3| - |B.6| ). Taking rri h = tnw, in Fig. 1 we 
compare the exact result (solid line) for the trace Tr (exp[—t/Cb os ] — exp[— 

(l.h.s. of eq. (|7.1|)) with several approximations (dashed and dotted lines), given by 
the first, the first two and all three terms of the r.h.s. of eq. ( |7.1| ). For low values of t 
we obtain excellent agreement between the numerical result and the approximations, 
as it should be. For large values of t the comparison does not provide a check of 
the numerical treatment. Here the approximations behave as some power law of 
t 1 / 2 while the exact result is dominated by the contribution of the negative and zero 
modes, given by + 6, which is also plotted in Fig. 1 (dashed line). 

Now that we have checked the reliability of the spectrum we can use it to calculate 
the desired quantities. These are the renormalized non-thermal contributions E FP 
and the temperature dependent parts Pbos a Fp(T) associated with the Bose-Einstein 
distribution factor. Both quantities have to be evaluated for the fluctuation operator 
/C b os and the Faddeev-Popov operator /C F p- 

To obtain the renormalized value of the non-thermal parts PbosFP we have to 
evaluate 77^ FP (A) in the limit of infinite proper time cutoff (A —> oo) (see eq. (|4. 1 1|) ). 
Numerically, however, we always have to work with a finite A to ensure the finiteness of 
the basis. For this reason the numerical parameters R (box size) and P max (maximum 
momentum) also have to be finite. In order to obtain the limit of infinite parameters 
we proceed as for the fermion non-thermal energy in [T7] : First we fix A and take 
R and P max large enough so that their further increasement would not change the 
result any more. This procedure is repeated with larger values of the cutoff A until 
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we can determine the limit E conv (A = oo) = E ren . We illustrate this method in 
Tabs. 1 and 2 at the example mn = rnw- The renormalization scale z/ ren is hxed 
according to eq. ( |4.25| ). With m t = 174 we obtain i/ ren = 2.02. In Tab. 1 we show 
results of E^ V (A) for hxed A = 4 and various values of R and P max . We find that 
for R — 12 and P max = 16 = 4 A the continuum limit P£°“ v (A = 4) = —6.28m w 
is reached with an accuracy of better than 0.2%. With the same method we obtain 
results for other values of A which are presented in Tab. 2. For large A the law 
r (A) = a + 6/A 2 is satisfied for the fit a = —5.95 mw and b = —5.7 mw- Thus 


v bos 


from the data in Tab. 2 we can extrapolate the result for infinite cutoff A and obtain 
T'bos = a = — 5.95mjy. Considering the possible error of the values for hxed A and 
the error of the extrapolation process we estimate an accuracy of better than 2% for 
the final result. For other Higgs masses the deviations can be slightly bigger but in 
general the numerical error for the non-thermal energy should be well below 5%. 

Table 3 shows results for the renormalized non-thermal energy for different Higgs 
masses including the contribution from the Faddeev-Popov operator. The renormal¬ 
ization scale z/ ren , hxed by eq. ( |4.25|) , is included in Tab. 3. As already mentioned in 
Sect. 4, for m H = 350 GeV eq. ([4.25| ) has no solution so that instead we choose z/ ren to 
minimize the difference between and the pole Higgs mass, which for this reason 
may deviate from 350 GeV by some value less than 10%. Moreover we give values for 
the classical sphaleron energy and for the fermion non-thermal energy. The latter was 
calculated for | heavy doublets with top quark masses between 150 and 200 GeV and 
9 + | massless doublets. We hnd that both the fermion and the boson non-thermal 
energy are significantly lower than the classical energy which is in accordance with 
the fact that, generally speaking, loop contributions are suppressed by a factor a 
relative to the tree contribution. Actually after the renormalization the non-thermal 
energy of the boson fluctuations about the sphaleron is small and negative while that 
of the fermions is larger and positive. However, we show below that the thermal part 
dominates the boson fluctuations and the sum of both parts has the same sign as the 
classical energy. 

The behavior of the non-thermal energies for low m# and large rn t can be described 
by simple scaling laws. For mt/mw > 1 the aggregate energy density of the Dirac 
sea is dominated by the square loop diagram in the external Higgs held and hence is 
proportional to V c (h$) 4 In (h$/z/ ren ) where h is the Yukawa coupling and $ the Higgs 
held of the sphaleron. To obtain the energy we have to integrate this value over the 
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space where the Higgs held differs from its vacuum expectation value, i.e. over the 
spread of the sphaleron. For small Higgs masses uin/mw 1 the size of the sphaleron 
fields roughly scales as m^ 1 since the asymtotic behavior for large radial distance r 
is dominated by the term e~( mH / mw ' >r . Therefore all spatial integrals and hence all 
matrix elements of the fluctuation matrices scale as mjj 3 . Thus for m# < my/ < m t 
the dependence of the zero temperature energies on m t and mu is roughly given by 


E r f en 

term 


ZTiren 

■^bos 


. zriren 
+ ^FP 


OC 

oc 


, AT mfln.(m t /v ien m w ) 

\ C q 

m H 

rr4y 

ni'fj 


(7.2) 

(7.3) 


These scaling laws explain the strong increase of the fluctuations for small and 
large m t , which is in correspondence to our numerical results. As mentioned above, 
we also found numerically that the boson and fermion non-thermal energies differ in 
sign. 

One finds a strong increase for small rri jj and large mt also for the thermal parts 
of the fluctuations, but here both fermionic and bosonic contributions are positive. 
For not too small temperature T the boson fluctuation energy is dominated by its 
positive thermal contribution so that the sum of the non-thermal and thermal parts 
is positive for both the boson and the fermion fluctuation. For small m# and large 
mt these sums are large so that they provide a strong suppression of the sphaleron 
transition rate. If rri h is small enough, this suppression prevents the erasure of the 
BAU. Thus, we see that the condition that the BAU should survive sets an upper 
limit on the Higgs mass. 

In order to obtain a quantitative result for this upper bound we still have to 

evaluate the thermal parts (see eqs. (B .41, 15.51)). To simplify notations we 

qmw - - 

will drop the subscript qmw in what follows. In principle we could evaluate them by 
a summation over the whole spectrum; numerically, however, it is preferable to sum 
only over eigenstates with low or medium energy and to use the expansion 


6{E) — <?oo + <?2 -jp + 


(7.4) 


(see eq. ( |B. 1| ) ) for the high energy part of the spectrum. To this end we take a 
smooth function F(E) with the properties F(E) = 1 for E < E a — Eb, F(E ) = 0 for 
E > E a + Eb and 0 < F(E) < 1 for E a — Eb < E < E a + Eb. Here E a and Eb are 
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fixed numerical energy cutoffs, usually we take E b = E a /2 so that E a is left as the 
only parameter. We calculate the thermal energies as follows: 


zrismall (n -?\ _ i 1 

- C/ bos,FPV 1 ) ~ 


1 _ p -Pqm w u> n 1 _ „-/3qm w ui° 

E FM In -E^n)ln 


/3qm w 
dE F(E) In - 


n 

D —(3qm w E 


/3qm w 


j3qmw 


q 2 


, 1 -F(E) , 1 
dE ——- In ■ 


^—PqmwE 


E 2 


(3qm w 


(7.5) 


The integrals in eq. (|7.5|) can easily be evaluated numerically, in the sum only states 
with u n < E a + E b appear. Now we have to check that numerically i7 sma11 is inde¬ 
pendent of the parameter E a . In Table 4 we show results of 3E^ aU (T) in the high 
temperature limit (see eq. (|B.15| )) for m# = rriw and several values of E a . We also 
give results for the contributions (the sum and the integrals of eq. (|7.5|) ) separately. 
We find that in the interval 3 < E a < 8 both the sum and the integrals in eq. 0 
drastically depend on E a but the result for (3E £™ a11 varies only by about 2%. For 
smaller values of E a the expansion eq. (|7.4|) is not good any more and for larger val¬ 
ues of E a the numerical accuracy of the spectrum decreases due to contributions of 
states with very large grand spin K. For other Higgs masses the variation of ^™ all (T) 
with E a can be about 5%. Hence for our calculations we choose 4 < E a < 6 and 
obtain E sman (T) with an accuracy of usually better than 5%. 

We are now going to compare our results to the ones obtained by Carson et. al. |I3 
and Baacke et. al. [TO] . There the expression 


In K = -PEl™ n (T) - (3E S ™ U (T) - 6 In 2 - In \u- 


small / 


(7.6) 


was evaluated in the high temperature limit T —> T c . In Fig. 2 we show the results of 
our work as well as those of [ 13 , !()] as a function of the Higgs mass m#. Our data 
are between those of |15|] and |16|, they agree with the ones of |1| up to 10%. Apart 
from numerical uncertainties one reason for the difference could be the renormalization 
scheme. We have performed the renormalization at zero temperature strictly, as it 
is usually done. This corresponds to a subtraction of the first term (£?oo) hi the high 
energy expansion (see eq. (|7.4| )) which is also the first term of the tadpole expansion 


3C]. In 116], however, all tadpole graphs except the term linear in T have been 


removed. The difference is then due to the higher order terms which are small for 
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high T but numerically not completely negligible. There is a larger deviation from 
the results of JT5|, only a qualitative agreement for the low m# behavior is found. 

Since the evaluation of the fluctuation determinants eqs. ( P . 2 6|^3728|) is a rather 
involved task one seeks for a good approximation procedure which is easy to handle. 
Before an exact calculation was performed, Carson and McLerran Jl3| applied an ap¬ 
proximation technique by Diakonov, Petrov, and Yung (DPY) [|14j to this problem. 
Later they calculated the boson fluctuation determinant exactly in the high temper¬ 
ature limit [b| and found that the exact and the approximative result deviate by 
several orders of magnitude. In this section we revisit the DPY method and explain 
how it can be used to get the the fluctuation determinants to a reasonable accuracy. 

Following [|L| we consider here the high temperature limit of the boson fluctuation 
determinant which can be written as an integral over spectral densities: 

InX bos = -f3 c E™ al \T c ) = -\J dE ( Q hos (E ) - ^ os ) ln(U 2 ). (7.7) 

where @ hos (E) is the spectral density of the boson operator, with the six zero and one 
negative mode subtracted. Using the identity 

roo dt / 
m a = / — e 

Jo t V 

we arrive at the proper-time representation: 

r°° dt 


—t 


—tot. 


) 


(7.8) 


In V., = 4 r dE (^(E) - r^(e-‘-e-“ a ) 

J 0 J 0 T/ 

= \ y | ( Tr exp[—f/Cbos] - 6 - exp(t|a;_| 2 ) - Tr exp) 

nbos r=r n 

(7.9) 


where Q^ s = 2a/y/n (see eq. (|B.8|) ). 


The idea of the DPY method is as follows [14]. The behavior of the integrand 


at small t can be established from the semiclassical expansion of the ’’heat kernel”, 
Tr exp (—tIC), see eqs. (|4.7| , |7T[) and (|B.3|) . Its behavior at large t is governed by 
negative and zero modes. Therefore, knowing the behavior of the integrand both at 
small and large t, one can approximate the integral of eq. ( |7.9| ) as a sum of small- 
and large-t contributions, separated by some parameter t 0 , 

rto 


lnXbos = / dt f{t) 


dt /low(f) T / dt /high(t) — fn kbosj 

/ 0 Jto 


(7.10) 
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where the separation parameter should be found from the requirement that the sum 
of the two terms in this equation is stable in to- Actually it means that to is a point 
where the small- and large-t approximations to the true integrand cross. If that does 
not happen the method fails. The more terms one knows from both sides, the better 
is the accuracy of the method. 

Using the heat kernel expansion for Tr exp(—t/C), we find the approximation for 
small t 

|4 ^ (7.11) 


f]ow(t) — -(&£ 2 +cf 2 — (7 + |u;_|Q + -(7 
At large t it behaves as 


u;_ 


/high(f) — 


1 (7e 


— at 


(7.12) 


2 V ( 

Knowing the coefficients a, b, c, and |u;_|, one can estimate the fluctuation determi¬ 
nant, eq. (O). using eq. ( |7.10|) . The result for three different values of the Higgs 
mass is presented in Table 5, together with the exact value of lnXbos- One observes 
that the accuracy is at the level of 10 to 15%, but that is a price one has to pay if 
one wishes to avoid a laborous computation of the exact spectrum. 

8 The upper bound for the Higgs mass 

The main issue of our work is the calculation of the sphaleron transition rate 7 = 
T jV according to the Langer-Affleck formula eqs. ( [3.11|. |5.12| ) including the classical 
Boltzmann factor, the fermionic and bosonic one-loop contributions and the Jacobian 
prefactors. We stress again that our calculation is not based on the high temperature 
limit but was done for arbitrary values of T. Therefore, we can compute the rate 
for the whole temperature range between zero and the critical temperature T c and 
perform the integration over T (see eq. ( |S .10 )) to obtain the ratio B 0 /Bt c ■ In Fig. 3 
we present the contributions —j3qE c i ass + In T (classical part), — /JU FERM (fermion 
loop), and —/3E B os (boson loop) of In 7 (according to eqs. (|5.12ft^TT^) ) for run = 66 
GeV and rri h = 100 GeV. It is convenient to take the parameter q = \J 1 — T 2 /T c 2 as 
independent variable rather than the temperature T itself. 

Qualitatively both pictures of Fig. 3 show the same behaviour, but we fold signif¬ 
icant quantitative differences. At low temperatures (large q) the main contribution 
to the classical part is the Boltzmann exponent ~/3qE c i asa which decreases with q 
roughly linearly. For T —> T c (q —> 0) the Jacobian prefactor In T ~ 71ng(T) —> — 00 



















dominates the classical part. Hence we find a maximum of the classical contribution 
to the transition rate at about q ~ 0.1. For large and medium q the suppression from 
the fermion loop contribution can be also rather large (in the case rn h = 66 GeV it 
becomes almost as large as the classical one) but it tends to zero in the high T limit 
q —>■ 0. The bosonic contribution is generally rather small and almost constant over 
the plotted range of temperatures. In the high temperature limit it does not disap¬ 
pear but it tends to some finite value. Therefore we see that in this limit which was 
assumed in previous works || §, [US]. [U]. |16[ the boson one-loop contribution is indeed 


the most important one while fermions decouple. Adding the loop contributions to 
the classical part, we obtain the total rate which has the same shape as the classical 
curve, in particular it also has a maximum. Hence if there was any significant baryon 
number violation after the electroweak phase transition, it must have happened in 
a short period around this maximum (remember that in Fig. 3 the logarithm In 7 is 
plotted while 7 itself enters the integral of eq. (|6.10|) ). We find that the piece of the 
curve for the total rate which is marked by a solid line contributes about 99% to the 
ratio log 10 (B 0 /Bt c ) which measures the washout of the BAU. 

Both the position of this washout area and the value of the maximum are strongly 
influenced by the loop corrections, especially we note that in this region the fermionic 
contribution which was neglected in previous works is quite essential. Below we inves¬ 
tigate the effect of the fermions quantitatively by computing the ratio B 0 /B Tr with 
and without fermion loop corrections and confirm the significance of the fermions. 

Comparing the two plots of Fig. 3, we find that the loop contributions are strong 
for low rriH so that in this case the rate is suppressed, while for large m# the fluc¬ 
tuations are rather weak. This is also documented in Fig. 4, where we have plotted 
the total transition rate for various m# and m t . In accordance to the scaling laws 
eqs. (|7.2| , |7.3|) we find a strong suppression of the transition rate 7 for small mn and 


large mt and a weak suppression for large rn h ■ This results finally for physically 
relevant m t in a small transition rate for small m# and a large transition rate for 
large m#. If the maximum value of 7 is small enough, the baryon number violating 
processes have happened so rarely that they have not affected the BAU. On the other 
hand a large transition rate means that the sphaleron transitions must have elimi¬ 
nated the baryon asymmetry. Thus, if we fix m t we obtain an upper bound for m# 
below which the asymmetry is conserved and above which we expect a dissipation of 
the baryon number. This is how we deduce our upper bound from the condition that 
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the BAU should survive the age of sphaleron transitions. 

Before we evaluate this conclusion quantitatively a comment on the limits of the 
approach is in order. Our calculation is based on the assumption that the Langer- 
Affleck formula is valid and our restriction to one-loop contributions is justified, 
i.e. that higher order corrections can be neglected. Both assumptions are reason¬ 
able for temperatures not too close to T c and if the fluctuations on the one-loop level 
are small compared to the classical part. The physics of the phase transition and in 
its direct vicinity, where perturbation expansion breaks down, is complicated and not 
well understood yet so it is difficult to decide at what temperature the framework of 
our calculation becomes inapplicable and what in this case could be a more adequate 
description. However, we can estimate the reliability of our model by checking if, first, 
in the washout region the loop contributions are not too big compared to the classical 
part, and second, if the onset of the sphaleron transitions, i.e. the left margin of the 
interval marked by a solid line in the corresponding curve of Fig. 4, is not too close 
to the critical temperature. Figs. 3 and 4 show that both conditions are well fulfilled 
for mu ~ 100 GeV but not so good for smaller rri h ; for run ~ 60 GeV the fluctuations 
are rather large so that the model on the one-loop level is probably not reliable. The 
conclusion of this restriction on the applicability of our technical framework will be 
drawn later. 

Furthermore, we assume that, before the transitions start, there is the same num¬ 
ber of baryons and leptons in the Universe, i.e. B — L = 0. In the standard model 
B — L is strictly conserved so that this condition will not change during the period 
of the transitions. If there were a primordial excess of e.g. antileptons, created by 
unknown forces which violate B — L before the electroweak phase transition, then the 
sphaleron transitions would increase rather than decrease the BAU. 

Let us also briefly comment on the connection between our critical temperature 
T c defined in eq. ( |5.9|) and the electroweak phase transition. As a consequence of the 
thermal renormalization the vacuum expectation value of $ becomes T dependent 
and vanishes for T —> T c . This looks like the behavior of fields at a second order 
phase transition. However, in order to obtain the true temperature and nature of the 
phase transition, one would have to include other terms, e.g. a term of the order of 
T<f> 3 into the potential. Since there is no consistent way how to perform calculations 
near the phase transition, where perturbation theory is not applicable, we decided to 
take only the numerically by far dominating term of the order T 2< L 2 explicitely into 
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the potential; other terms, like the T $ 3 one, are considered in the quantum correction 
E smaR . Thus our critical temperature should be seen as a mere definition which needs 
not necessarily coincide with the temperature of the true transition, neither do we 
imply that the phase transition is of second order. 

Knowing the transition rate 7 as a function of q for a fixed Higgs mass m# it 
is now possible to perform the integration in eq. (|6.10|) numerically. The result for 
the ratio log m (B 0 /B Tc ) is plotted in Fig. 5 for m t =150, 174, and 200 GeV. For 
comparison, we also performed the calculation without considering fermions. In this 
case we did not use eq. (|4.25|) to fix the renormalization scale, but a corresponding 
equation which follows from the Higgs propagator with boson fluctuations included 
instead of fermions. Here we obtain z/ ren ~ 1. 

All the curves start at zero for small Higgs masses which means that the BAU 
survives completely. If we increase m#, the fluctuations become weaker so that the 
transition rate increases. Hence the ratio Bq/Bt c suddenly begins to fall, and within 
a short interval it drops by 20 orders of magnitude. The bigger the top quark mass is, 
the larger are the fluctuations, and hence the region where this decrease takes place 
is shifted to larger Higgs masses. For a Higgs mass beyond this region the survival 
of the BAU is ruled out, irrespectively of its initial value Bt c immediately after the 
phase transition. Let us assume that this initial value is such that in order to explain 
the present day BAU we have to demand B 0 /B Tc > 10 “ 5 [Q (we see, however, from 
Fig. 5 that this value is not important since a change by many orders of magnitude 
alters the upper bound by only a few GeV). For m t = 150 GeV we obtain m# < 60 
GeV, for rn t = 174 GeV the upper bound is at 65 GeV, and for m t = 200 GeV the 
BAU survives if < 71 GeV. At any rate, for all physical choices of the parameters 
m t and B 0 /Bt c the upper bound for the Higgs mass is found in the range between 60 
to 75 GeV. 

The calculation without fermions leads to a qualitatively similar picture, but the 
erasure of the asymmetry happens already at much lower Higgs masses. Assuming 
Bq/B Tc > 10 -5 we would obtain an upper bound for rn h of only 49 GeV being close 
to the bound found previously by Bochkarev and Shaposhnikov !, !• The large 
difference between this value and the bound of 65 GeV which we obtain with fermion 
fluctuations for rn t = 174 GeV again confirms their significance. 
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9 Summary and Conclusions 


The present paper investigates the fate of the baryon number asymmetry in the 
Universe (BAU) after the electroweak phase transition (T < T c ). It is assumed that 
the asymmetry as such originates from baryogenic processes before or during the phase 
transition with a net result of B + L ^ 0 and B — L = 0. It is furthermore assumed 
that in the broken phase (T < T c ) the Minimal Standard Model with one Higgs 
doublet holds. Since the Standard Model does not conserve the baryon number due 
to possible sphaleron transitions, today’s existence of an asymmetry of about 1CU 10 
baryons per photon implies certain dynamic conditions right after the phase transition 
which prevent too fast a “wash-out” of the baryon number [y]. In the present paper 
the baryon number transition rate is evaluated in the one-loop approximation around 
the classical sphaleron solution (hedgehog). The higher loop effects are partially taken 
care of by an exact treatment of the “Debye mass” terms, ~ $ 2 T 2 . It is assumed that 
the Langer-Affleck formula holds and that no baryons are generated in the broken 
phase. 

For all temperatures below the critical temperature T c the one-loop calculations 
are performed numerically in the limit of vanishing Weinberg angle. In fact, the 
baryon number transition rate depends on the classical sphaleron energy, the deter¬ 
minant of the fermionic fluctuations, the determinant of the non-zero bosonic fluc¬ 
tuations, the energy of the negative mode and the normalization factors of the zero 
modes. While the sphaleron energy and the zero and negative bosonic modes have 
been calculated previously in the literature [||, the evaluation of non-zero bosonic 
modes has been performed only in the high temperature limit with somewhat contro¬ 
versial results |T5|, |T(|. In this context the present paper shows the first calculation 
of the boson determinants for finite temperatures (the fermion determinant at arbi¬ 
trary temperatures was previously computed by the same authors JlTj). It turns out 
that all above contributions to the transition rate are more or less equally important 
and must be evaluated at finite temperatures in order to obtain, within the given 
conceptual frame, an accurate calculation j[U] . 


The actual numbers basically depend only on one unknown parameter, namely 
the mass of the Higgs boson m#. In fact, the dependence of the baryon number 
transition rate on the Higgs mass is extremely strong. Both bosonic and fermionic 
fluctuations above the sphaleron barrier help to preserve the baryon asymmetry in 
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the Universe. They prevent a fast erasure of the baryon excess provided the mass 
of the Higgs boson is less than some upper bound, while for larger Higgs masses 
the sphaleron transition rate becomes large and the asymmetry would be eliminated. 
The value of this upper bound depends on the mass of the top quark, ranging from 
about 60 GeV for m t = 150 GeV to 71 GeV for rn t = 200 GeV. These results are 
obtained in the Minimal Standard Model with only one Higgs doublet. They assume 
a theoretical frame characterized by the applicability of the Langer-Affleck formula 
and the restriction to one-loop calculations, with a partial resummation of higher 
orders. These assumptions are only justified if the baryon number violating processes 
do not happen immediately after the electroweak phase transition where the loop 
expansion breaks down. This means the position of the maximum of the transition 
rate 7 should be not too close to the critical temperature. Moreover the quantum 
loop contributions at the maximum should be small compared to the classical terms. 
We find that both conditions are well fulfilled for Higgs masses ran ~ 100 GeV while 
this is not the case for small Higgs masses below about 60 GeV. Those Higgs masses 
are, however, ruled out by experiment [|3T], [£|. Thus we arrive at the following 
conclusion: If the Higgs mass is in the range between about 60 and 100 GeV, the 
Minimal Standard Model could be able to account for the survival of the BAU, either 
within the formal framework we used and a suitable top quark mass or by effects 
outside this formalism, e.g. higher loop contributions. If it is found above 100 GeV, 
there is only a little chance to explain the present BAU within the MSM since in this 
case the application of our framework is rather safe and predicts the complete erasure 
of the BAU. A possible escape could be an extended model with two Higgs doublets, 
following from supersymmetric models. 
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Appendix A 

In this appendix we describe how we solve numerically the eigenvalue problems 


/C hns T hns = u/'F 


bos 


/Cfp'&fp — ctv'b 


FP 


(A.l) 


for the boson fluctuation and the Faddeev-Popov operators. We construct a finite 
basis for the fluctuations in which the operators can numerically be diagonalized. 
Partially this technique has been developed in the context of the chiral quark model 
and employed for the diagonalization of the fermionic fluctuation matrix 


17], which ensures consistency between the calculations of bosonic and fermionic loop 


corrections. 

The fluctuation vector Tbos consists of nine gauge field components af (a, i = 
1,..., 3) and four Higgs field components ip M (/i = 1,..., 4) while Tpp contains only 
three components which we denote by Qq. Hence the eigenvalue equations read: 





= u 



T ab n b — } 2 n a 

J lIq — UJ CLq , 


(A.2) 


where the matrix elements of the fluctuation matrix IC^ os are given by (see eq. ( 3.21 )) 
9$ = 5 ij 5 ab (-d 2 + A c k A c k + \$ fl $ fl )+2£ abc Fc 


v 


= 


+ S,j {e^(d k AD + 2e abc A k d k - A‘ k A b k ) , 

V {-9 2 + l A ‘ A ° + IV, + - 4 )) 

+ + \{v\ - , 

= ls abc rj^%A b ~ , 

and the Faddeev-Popov operator by (see eq. (|3.14l) ) 

T ab = 5 ab (_ d 2 + A c A c + J + £^(3^) + 2 S^Afdi ~ A“A b 


W c 

' Y n 


(A.3) 


(A.4) 


For the static classical sphaleron fields we assume the spherically symmetric hedgehog 
ansatz: 


Aaf \ 1 — A(r) N-S(r) C(r) 

(rj ^aijUj " 1“ \Vai " 1“ 


<3h(r) = 2 mG(r), $ 4 (r) = 2 H{r) 


(A.5) 


with the given five radial functions A(r), B(r ), C(r), H{r ), G(r). In principle one of 
these functions can be eliminated by a gauge transformation. Although this would 
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lead to a significant simplification of the numerics we will not perform this step but 
rather stick to the general gauge with five functions since our numerical procedure 
only works if the classical fields are continuous functions at zero and infinity, i.e. that 
they take their vacuum values there. This is only possible for a non-vanishing C 
held, which increases the numerical effort, but on the other hand allows to check the 
invariance of all quantities under gauge transformations. 

To exploit the spherical symmetry and to construct a finite basis in which /Cbos 
and /Cfp can numerically be diagonalized we consider a four (= 3 + 1) dimensional 
reducible SU(2) representation with generators: 


Si = 


< 0 0 0 o \ 

0 0 -i 0 
0 i 0 0 

0 0 0 0 


, S 2 = 


( 0 0 i 0 \ 

0 0 0 0 
-i 0 0 0 
0 0 0 0 


, s 3 = 


/ 0 -i 0 0 \ 
i 0 0 0 


0 0 0 0 
0 0 0 0 


(A.6) 


A basis of (“spin”-) eigenstates \S S 3 ) of S 2 = Sf + S% + S$ and S 3 is given by 

( o\ ( ~i\ f Q \ fi\ 


|0 0 ) = 


0 
0 

V 1 / 


’ |11) = 72 


1 
0 

V 0 I 


, | 10 ) = 


0 

i 

v°y 


’ 11 ~ 1)= V2 


1 
0 

v°y 


(A.7) 


The indices i,j will always refer to coordinates of these spin eigenstates, for example 
|1 l)i=i = — i/ a/2, |0 0 )j=4 = 1. Similarly we define a four dimensional SU (2) “isospin” 
representation; the operators 7\, T 2 and T 3 and the eigenstates |TT 3 ) look exactly as 
the corresponding ones of the spin representation. Here the coordinates are referred 
to by the indices a,b and Moreover we use the basis \LL 3 ) of the angular 

momentum operator to describe the spherical space dependence of the fluctuations, 
with the property 

<fi| LL 3 )=i l Y LL 3 (SI). (A.8) 


The “grand-spin” operator defined byK = J + T = L + S + T commutes with the 
fluctuation operators of eq. ( |A.1| ) . Therefore eigenstates of K 2 and K :i form a proper 
basis for the diagonalization procedure. We couple the eigenstates \L L 3 ), \S S 3 ) and 
| TT 3 ) to eigenstates of K 2 and Ji 3 : 


I K, K 3 ] T, J, S, L)t = E Cf £ TT3 CHlss, \SS3)r | TT 3 ) a \LL 3 ). (A.9) 

L 3 ,S3,J3,T 3 
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For the vacuum fluctuation operator, i.e. in the case of no external field, we can 
solve the eigenvalue problem analytically; the dependence of the fluctuations on the 
radial coordinate r is in this case given by spherical Bessel functions. We take these 
solutions as the basis for a numerical diagonalization of /C in the non-vacuum case. 
To this end we define states | p\ K, K :i ; T, J, S, L ) by 

(r | p- K, K,- T, J, S, L) = A f j L (pr) <fi| K, K :i] T, J, S, L ). (A.10) 


Here the momentum p is a continuous variable, and J\f is a normalization factor 
specified below. In order to get a finite basis we have to discretize the momentum 
and to restrict its allowed values to a finite number. With large enough numerical 
box parameters R and P max we demand: 


where 


Ji(PnR') = 0 , 

I = I{K, J, S) = 


P a < Pn 


(A.ll) 
(A-12) 


K for S' = 0 
J for S' = 1 

In the case S = T = 1 we have three discretization conditions for fixed grand-spin K 
instead of one, yielding three sets of momenta p^ +1 , p%, and p^ _1 . This extension 


of the usual construction j34| is necessary to ensure the orthogonality of the basis 
states; it has already been used and checked in [^] . We obtain 

r R 


K.Kn I dr r2 JriPn r ) 31 (Pm r 
Jo 

rR 


= KiK \ dr r 2 Jl±l(Pn r ) jl±l {pL r ) = s nm , (A. 13) 

Jo 


if the normalization factor is chosen as 


2 . ' T —' — i 


^n = \h^ Jl±MR) 


R 3 


(A. 14) 


hence our states are orthonormal: 


(p„; K, K 3 - T, J, 5, L \p.P K\ K’p r, ./', S', U) 


3nm d KK’ d K 3 K' 3 d TT' d JJ' d SS' d LL' 


(A.15) 


For fixed values of K = 0,1, 2,... and K 3 = —K ,..., +K we can write down the 
following set of basis states for the fluctuations: 




1 , a I 

bos 


r = 


' <( r ) ^ 


f (r\P J n ,K, K 3] 1, J, 1, L)^ 



Pa{ r) 

= 

0 

= 

0 

V P4(r) j 


l 0 J 


l 0 ) 
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for J = K — 1, K , K + 1; L — J — 1, J, J + 1; n = 1,..., N( J); 



' <(r) N 


( 0 ^ 


( 0 \ 

'Fbos(r) = 

7 >a(r) 

= 

( r Pn ; K s; 1, l, 0, L)% 

= 

<rfc:>4 a 


^ ^4(r) ) 


{ 0 ) 


l 0 J 


for L = K - 1, K, K + l] n = 1,..., #(#); 


1 <(r) N 


f ° ^ 


{ ° 1 

^o(r) 

= 

0 

= 

0 

V M r ) J 


V < r Pni K, K 3l 0 ,K, 0, K)\ ) 


v ( r l^bos)4 ) 


for n = 1,..., N(K)] 


^pp(r) = ag(r) = (r |p*; K, K 3] 1, L, 0, L)% = (r |T“ P >“ 

for L — K— 1, K, K+ 1; n = 1,..., 7V(iT); (A.16) 


where iV(/) is the number of allowed momenta p T n , see eq. ( A.11 ). The index a enumer¬ 
ates the basis states of the three groups for the fluctuations and of the Faddeev-Popov 
matrix. For K = 0,1 not all of these basis states exist. The total number of states 
for T bos is given by 3 N(K + 1) + 7N(K) + 3N(K — 1) for fixed K > 1 and fixed 
K :h for K = 1 it is 3N(2) + 7iV(l) + iV(0) and 3iV(l) + 2N(0) for K = 0. For the 
Faddeev-Popov matrix we have 3N(K) basis vectors for K > 0 and A r (0) for K = 0. 

We show below that due to the spherical symmetry of the sphaleron the operator 
K , bos is block diagonal in K and K 3 , i.e. basis states with different K or K 3 do not 
mix. Moreover the blocks for different K 3 and the same K are identical, so that for 
each K only one matrix has to be diagonalized, and the resulting eigenvalues are 
(2A'+l)-fold degenerate. The dimension of this matrix is given by the number of the 
above basis states. The same holds for the Faddeev-Popov operator /Cpp- 

The remaining task is to calculate the matrix elements of the operators in the 


basis (|A.16|); i.e. if | v h b os Q1 ) and | v h^ Q2 ) are basis states given by eq. (|A.16|), we need 


• A2,02 


to know the element (Tb^ Ql |AA os |Tb 2 g“ 2 ). F° r this purpose we have to express the 
matrices in terms of spherical tensor operators so that the spherical part of the matrix 
elements can be evaluated analytically. Apart from S and T given in eq. (|A.6|) we 
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need operators P + and P , acting in spin space, which we define as 


P+ = 


( 0 0 0 1 \ 
0 0 0 0 
0 0 0 0 
10 0 0 


and 


pr = 


(0 0 0 
0 0 0 
0 0 0 
Y % o o 


-i\ 
0 
0 

0 / 


p+ — 
- r 2 — 


( 0 0 0 0 \ 
0 0 0 1 
0 0 0 0 
0 10 0 


P-2 = 


(0 0 0 
0 0 0 
0 0 0 
0 i 0 


0 

0 


Pt = 


i 0 0 0 0 \ 
0 0 0 0 
0 0 0 1 
0 0 10 


P7 = 


( 0 0 0 0 \ 
0 0 0 0 
0 0 0 -i 
0 0 i 0 


(A.17) 


(A.18) 


It can easily be checked that like S these operators are spherical vector operators, 
i.e. [Ki , Pp] = i£ijkPk ■ Moreover we define two scalar operators by 


I,= 


( 1 0 0 0 \ 
0 10 0 
0 0 10 
0 0 0 0 


1 0 0 0 0 N 

0 0 0 0 

0 0 0 0 

v 0 0 0 1 J 


(A.19) 


In the isospin space we need the corresponding operators which we denote by and 
It, *t- 

Using the relations 


hj — ( Is)ij ) fi ab — ( 7 t )“ 6 , Spu — (It + , 

Sijk = e abc = i(T c y b , v% = i(T a + QZU, 

A a k A\ = (A c k A c k I T - A c k A d k T d T c ) ab , W r a b = (W k rk I T - W r ^T fc )“ & , 

— (&k$klT + — $l$kTlTk + )nv , 

(Pr )*4 (Pr ju (^P'r )i4 Pr )4 i ^ri , 

(Qtr = (Qtr = w~ k y 4 =-w~ k ) 4a = s ka ( a . 20 ) 


we can rewrite the fluctuation matrices: 

S$ = {W-d 2 + + IsT c T d AlAt 

^ ab 

+ i T c (2 I s A c k d k + I s (d k A c k ) + i e klm F c kl S m ) , 

) ij 
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= {i s (/T + ir) (-9 ! +K”A‘+ i*r% + yj*,*r-±)) 

+ i S (T a + Q:)i (A“9i + !(M?)) 

+ is 3 f s'i — 1) + *&4$4*T — Ql^kTlTk + Qk^iQt )| i 

w“ = {i(w; t / T -w5r 1 T i ) (j^ + iF-) 

a 

+ 3 W rfWt +iQt)(Pr + +iPr)\ , 

J i4 

M*. = {i(w; t /r-W5TiT 1 ) (P+-iP-) 

+ iwUQt-iQ*)(p?-ipn\ , 

T ab = {isIr(-d 2 + ^^^ + isT c T d A^ 

\ ab 

+ i s T c i(2A c i d i + {d i A?))\ . (A.21) 

J 44 

Now we can find matrices JC^ os and K-pp with the property 


<* 


Ai,ai 

bos 


(^FP 


| A^bos 

ai |/C 


Kr a > = 




FP 


1 ^) = 


(^ p \JC FP \ni) 


(for the definition of the states | \ 4 os) and | d'pp) see eq. ( |A . 1 6|) ): 


(A. 22 ) 


/Cbos = I T I s (-d 2 + ^^,) + I s T c T d AlA d k + I s T c {2A c k id k +t(d k A c k )) 

- e klm F c kl S m T c + i s (I T + i T ) (-a 2 + |A“A“) 

+ isIr(Q + §4) $ 4$4 + |4$k$* -14) 

+ + 14) + 14^4^4 -14) 

+ + Q a ) + |(f 9jA“)) 

+ i(4 - 1) (4Q^$4 - isT,r fc $ fc *,) + J T P+W r fc fc 

- i £Wm W r fe ; r m p- - i (w* + W' fc ) TiT k P+ 

+ \^rA {QtPr + ~ Qk P r) , 

Kfp = i s /T(-a 2 + i$ /i $ M ) + i5TcT d A l c Af + i s r c (2^i9 i + i(a i ^)). 

(A.23) 


Onr next step is to plug in the hedgehog ansatz (|A.5|) so that the matrices /Cbos and 
^fp can be expressed through the profile functions A, B , C, H , G, and the spherical 
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vector and scalar operators. With the hedgehog ansatz and 


&i f^i q t r ^ijk^jLk (A.24) 

we obtain after a long and tedious calculation: 

ic bos = (-ft-f £ + % + G 2 + H 2 + ± ((l- A) 2 + B 2 )y T I s 

+ (G 2 - B 2 - (1 - A) 2 ) I s (n • T) 2 + ^ (1 - A) I s (T • L) 

+ 4, (2 rC£ + rC' + c) I s (n • T) + f I s (T • (n x L) - i (n ■ T)) 

+ £(l - A 2 - B 2 + rA' + BC) (n • S) (n • T) 

+ £ (rtf' - AC) n-(SxT)-^ (rA' + BC) (S ■ T) 

+ (-^3 - f ^: + + 2^2 ((1 - A) 2 + 5 2 + ^-)^s(ir + *r) 

+ (V 2 + \v 2 h {h 2 - i) + y H G 2 y s i T 

+ ^G 2 + |^ 2 (G 2 — 1) + |+ (1 — ^i) G 2 i$ (n ■ T) 2 
A(l-A)i 5 (T + Q-)-L + A 5 (2rG| ; + rG' + G) i s n-(T + Q-) 

+ £^s ((T + Q-) • (n x L) - i n • (T + Q~)) 

- (1 - u 2 h ) HG i s (n ■ Q + ) + i (GG + 2r H') I T (n ■ P + ) 

- i (G + GA - tf£) (T • P") + i (if - ifA - £G) n • (T x P~) 

+ ± (G + GA - HB + HC - 2rG') (n ■ T) (n ■ P") 

+ ±{H - HA- BG ) n (Q + x P+ - Q x P ) 

- A (G + GA - HB) (Q + P + - Q P-) 

+ ±(G + GA-HB + HC - 2rG') 

((n • Q + ) (n ■ P + ) - (n • Q“) (n • P~)) , 

Afp = (-ft-f i + ^ + G 2 + H 2 + ^( y (l-A) 2 + B 2 )y s lT 

+ A (G 2 - B 2 - (1 - A) 2 ) i s (n • T) 2 + £ (1 — A) (T • L) 

+ A ( 2 rG|: + rG' + G) *s (n ■ T) + f (t • (n x L) - i (n • T)) , 

(A.25) 


where we dropped the argument r of the profile functions. Now it is easy to see that 
all the spherical operators which turn up in eq. (A.25) are scalar operators in the 
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sense that they commute with the grand-spin K 2 and K 3 so that 


/C bos , K 2 

= 

K-'bos 5 -^-3 

= 0, 

/Cfp , K 2 

= 

Afp , A3 


= 0 . 


(A.26) 


This is the reason why the matrices can be diagonalized in each K sector separately, 
as mentioned above. 

For the numerical diagonalization one has to evaluate these matrices in the basis 
(|A.16| ). The radial part of a matrix element leads to a numerical computation of a 
one dimensional integral. The angular part, however, can be evaluated analytically. 
The most direct and easiest way to do this is to employ the Wigner-Eckart theorem 
and perform the summations over the Clebsch-Gordan coefficients with the help of 
a program like Mathematica. Instead of writing down the complete result of this 
angular integration we demonstrate the procedure with an example. We consider the 
operator I (GC + 2 rH') It (n • P + ) which is part of Ab os and calculate the matrix 
element 

(p J n, K, K 3] 1, J, 1, L | i (GC + 2 rH') I T (n ■ P+) | K, K 3 , 1, I/, 0, L') 

= [ dr r 2 j L ( p J n r) ± (GC + 2 rH') j L , (p*) 

Jo 

K, K 3] 1, ,/, 1, L I T (n • P + ) | K, K 3 - 1, I/, 0, V ) . (A.27) 

The integral over the radial coordinate r has to be evaluated numerically. The angular 
matrix element is independent of K 3 so that we can put K 3 to zero. We obtain: 

( K, 0; 1, J, 1, L I T (n • P + ) | K, 0; 1, L', 0, L' ) 

i 

E \ ' \ ' ( -i \mr-iK 0 r~iJ J3 0 

1 -v U J j 3 ,it 3 ( -'ll 3 ,is 3 

L3,S3,J3,T3 L' 3 ,T ,^ m =—1 

•(1T 3 |/ t |1T') (15 3 |P ( +)|00) (L L 3 \ n ( _ m) | L' L' 3 ). (A.28) 

Here the operators PW and ri(_ m ) are spherical, not cartesian components of P + 
and n. The matrix elements in the last line can subsequently be evaluated using the 
Wigner-Eckart theorem, e.g. 


(ISsI Pf m) 100) = ( 1 IIP + II0). 


V3 


(A.29) 


Hence for the calculation of the matrix element in eq. (|A.27| ) we need to perform a 
numerical integration, a summation over Clebsch-Gordan coefficients, and to know 
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the following reduced matrix elements: 


(L 11 L 111/) 
(L 11 n 111/) 

( L 11 n x L 11 L') 


(S' 11 S 11 S') 
<S||/ s ||S') 
(5'HisllS") 
( S 11 P + 11 S') 
(S 11 p-11 S') 


hu \fL(L + 1)(2L + 1), 

8\L—L'\, 1 ( — *) \J\(L + L' + 1) , 

(-(L-1)VL lor // - L- \ 

\(L + 2)yfL + l forL'= L + 1 , 
10 otherwise 

8 S s' \jS(S + 1)(2S + 1), 

a/ 3 551 <5g/i , 

<5so ^S'O , 

—i x/3 (5si <Vo + ^so &s'i) j 

—x/3 (5gi 5g'o — ^50 <5g'i) • 


(A.30) 


For the correspoding reduced matrix elements of the isospin states one simply has to 
replace S by T and P by Q. 


Appendix B 


We discuss here the spectral densities 6 (E) defined in eq. ( |5.1|) , especially their asymp¬ 
totic behavior for large E. It is easy to see that at large E one can expand the spectral 
density in a series 

OO 

6 (E) ~ 6oo + £ 6 2 n E~ 2n . (B.l) 

n= 1 

In order to calculate the values of the coefficients 

lim 6 (E) = 600 and lim (6(E) — Q^E 2 = 62 (B.2) 

E—>oo E —>oo V / 

we use the small t expansion of 

F(t) = Tr (exp[—f/C] - exp[-t/C (0) ]) = ]T dE 6 (E) e ~ tE2 

discrete ^ 

levels 

= at ~ 1/2 + bt 1/2 + ct 3/2 + ■■■ . (B.3) 


The corresponding coefficients can be easily calculated using gradient expansion, fol¬ 
lowing the renormalization procedure eqs. (|4 . 7| — 14 .1 OP one can read off 


app — 


32tt3/2 


J ^($*$-4) , 
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^bos 
^ferm 
bpp = 
bhos 

^ferm 


327r 3 / 2 


(4 + u 2 h ) J d 3 r (4> f <I> - 4) , 


-^Nrfjd 3 r(#'i>-4) (B.4) 

TS+ / <i 3 r (-i(i4 “) 2 + £(*'* - 4 ) 2 + |(#'4 - 4)) , 

Tab-/(1r(+) 2 + lit 4 - 3<4 + bKi’i - + , 

+ |(8 - 3 vl + b)(4>*4 - 4)) , 

T5+r / dh (j N„(N C + 1)(F “) 2 + I NrfW - bXi 2 * - 4 ) 2 , 

+ \Nrf(i V 3 - V l)(m-i)) , (B.5) 

where we made use of eq. Q3.8|) to eliminate the term (A^A(A5*) in &bos and frf erm . 


In Section 7 we also need ||13 

E 3 r 


Cbos 


1 


384 tt 3 / 2 


o 2 ttt ci T?a | 28 _ abcjpa rpb jpc 

+ 15 £ jk^U 


+ \ (-34 + 93) ($t|>) /•)'/- + | (54 - 44 + *f) 4($t$)] 2 

+ i (4 + 284 + f) 4^) (A*) f (A*) 

+ + (154 + 214 + 184 + 48) 4 f <h - 4) 3 
+ | (274 + 574 + 364 + 144) ($1$ - 4) 2 


+ 9(4 + 24 + 4 + 8 ) ($*$ - 4) . 

Using (|B.3| ) one can show that 

lim Vi F(t) = 14^ , 

which implies 

Poo = 2a/va ■ 

To calculate £+ let us introduce the function R(E ) such that 

R\E ) = Q(E) and i?(0) = n D = number of discrete levels , 


(B. 6 ) 

(B.7) 

(B. 8 ) 

(B-9) 


n 


FP 


= 0, < os = 7, 4 rm = 1. Using partial integration, we can write 


Vi F{t ) — a _ 1 , 
t “ Vt ' ^ 


—tuj* 


discrete 

levels 


+ J™dE (6{E)- 6 00 ) e~ tE2 


1 

Vt 

1 

Vt 


E 

discrete 

levels 


—tbJ J 


+ (R{E) 


Q*>E\ e 


„~tE 2 


+ 2 Vi f™ dE (. R(E ) - Q^E) E e ~ tE2 


E < 

discrete 

levels 


—tcu 


1 ) + 2 Vi J™ dE (R(E) - Q^E) E e ~ tE2 


(B.10) 
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The first term vanishes with t —» 0, and for the second one we can use the same way 
as above with (R(E) — PooE^J E instead of 6 (E) and find that 


Jim (R(E) — QoqE^) E 



Using l’Hospital’s rule we finally find 

6 2 = Jim ( 6 (E) - 6 oo) E 2 



(B.ll) 


(B.12) 


From eq. ( |B . 1 1| ) we can deduce another interesting result; Jim (^R(E) — 6 oc E^j = 0 
yields 


J dE (j 6 (E) — £?oo) = Jim (^R(E) — 6 ooE — R( 0 )) = ~R( 0) = —n D . (B.13) 

Finally we show that the contributions i ? sma11 to the transition rate 7 are finite 
in the high temperature limit, i.e. q —> 0. In the fermionic case we obtain by partial 
integration 


^small 

7 ferm 


qmw 


ro o , v 

-i / dE (Q ierm (E) - 6 f £ m ) ln(l + e - fiqmwE ) - J < rm In2 

J 0 

/ \ r°° B ierm (E) 

— „ferm\ i„ o \ J 


( 0 ) - < m .) 


In 2 — qmw / dE ■ 


6 f S m E 


X _|_ e f3qm w E 


q ^0 


0 . 


(B.14) 


For El ™? 11 and UjpS a11 we obtain 


^bos a,AAU - -^FP 
^gismall 


(3qm w E 


qm w 


= ± 




ro o , , 1 

l3qmw 

J°° dE (e- (E) - 6^) In E (finite) (B.15) 


The last line holds since this is true for the integral J E with arbitrary upper bound E, 
and the rest /^° vanishes with growing E due to the behavior of 6 (E) — 600 , eq. (|B. 12 |). 
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R 

10 

10 

10 

12 

14 

P 

1 max 

12 

14 

16 

16 

16 

S£“ v (A = 4)/m„, 

-6.25 

-6.29 

-6.29 

-6.28 

-6.28 


Tab. 1: P{;°“ v (A = 4) for different values of the numerical parameters R and P max . 
The result shows that R = 12 and P max = 16 = 4 A are large enough to ensure that 
the continuum limit is reached. The mass parameters are = mw = 83 GeV, m t = 
174 GeV, z/ ren = 2.02. 


A 

2 

3 

4 

4.5 

5 

5.5 

6 

£K V (A )l™w 

-6.85 

-6.47 

-6.28 

-6.22 

-6.18 

-6.14 

-6.11 


Tab. 2: E^ V (A) for various values of A. Using the behaviour P£°“ v (A) = E^+b/A 2 
for large A, we find by extrapolation the continuum limit = — 5.95 mw- The mass 
parameters are m# = mw = 83 GeV, m t = 174 GeV, z/ ren = 2.02. 


rriH [GeV] 

50 

66 

66 

66 

83 

100 

150 

350 

m t [GeV] 

174 

150 

174 

200 

174 

174 

174 

174 

^ren 

2.07 

1.74 

2.05 

2.36 

2.02 

1.98 

1.76 

2.11 

E bos/ m W 

-9.74 

-5.40 

-7.22 

-9.32 

-5.95 

-5.09 

-3.64 

-5.69 

Epp/rriw 

1.36 

0.48 

0.75 

1.10 

0.46 

0.29 

0.08 

-0.01 

(E’C + E?$)/m w 

-8.38 

-4.92 

-6.47 

-8.22 

-5.49 

-4.80 

-3.56 

-5.69 

Efojmw 

26.94 

10.35 

16.91 

26.81 

12.11 

9.40 

5.87 

4.05 

F/class/ TTlW 

96.94 

99.60 

99.60 

99.60 

101.94 

104.08 

109.27 

121.67 


Tab. 3: The renormalized non-thermal energy of the boson fluctuations E of the 
Faddeev-Popov operator E§ p, and the sum of both for various values of the Higgs 
and the top quark mass. For comparison the classical sphaleron energy E class and 
the fermionic non-thermal energy E\^} m are included. The renormalization scale is 


determined according to eq. (|4.25|). 
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E a 

2.0 

3.0 

4.0 

6.0 

8.0 

E b 

1.0 

1.5 

2.0 

3.0 

4.0 

1st line of eq. (|7.5|) (sum) 

2.54 

10.95 

21.93 

49.95 

83.95 

2nd and 3rd lines of eq. (|7.5|) (integrals) 

3.64 

-4.20 

-15.07 

-43.10 

-77.13 

ft^r"(r c ) 

6.18 

6.74 

6.85 

6.85 

6.82 


Tab. 4: /9 c ^“ all (T c ) and its contributions for several values of the numerical param¬ 
eters E a and E b . The contributions strongly depend on E a and E b , but Eb™ a11 is ver y 
stable in the range 3 < E a < 8. 


run [GeV] 

66 

83 

125 

In /Y.bos 

-11.66 

-6.85 

-1.96 

In Xbos 

-12.91 

-7.66 

-2.27 


Tab. 5: Exact and approximate results for the high temperature limit of the boson 
fluctuation determinant for various m The exact values lnXbos are determined 
by a summation over the spectrum of eigenvalues of the fluctuation operator, the 
approximate values lnXbos are obtained with the DPY-method ||T4 
accuracy of about 10 to 15%. 


One finds an 
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Heat kernel approximations 



t 


Fig. 1 : Exact (solid line) and approximate values (dashed and dotted lines) of the 
heat kernel Tr (exp[—t/Cbos] — exp[—in dependence of the proper time param¬ 
eter t. The exact result is obtained with the discretized spectrum of the fluctua¬ 
tion operator Afbos, the approximations are the first three orders of the expansion in 
eq. (7.1). For low t we obtain excellent agreement. The large t behaviour of the heat 
kernel is governed by the negative and zero modes (dashed line). 




Boson determinant in high T limit 





Fig. 2: In k = — /h^b™ all (rc) — ^cEpp^iTc) — 6 In 2 — In |cu_ | as a function of the Higgs 
mass. We compare our results (solid line) with those of Baacke et. al. [16] (dashed 
line) and Carson et. al. [15] (dotted line). 




Contributions to transition rate 



0.00 0.10 0.20 0.30 0.40 0.50 


q 



0.00 0.10 0.20 0.30 0.40 0.50 


q 

Fig. 3 : Contributions to the transition rate In 7 for fixed mu = 66 GeV, m t = 
174 GeV and ra# = 100 GeV, m t = 174 GeV as a function of the parameter q = 
yjl — T 2 /T c 2 . The total rate has a maximum close to the critical temperature. The 
region which mainly contributes to the integral eq. ( 6 . 10 ) is marked by a solid line. 







The sphaleron transition rate 

1.0 0.999 0.99 T ! Tc 0.95 0.9 



q 


Fig. 4 : The sphaleron transition rate In (cil m w) for various values of mu (given with¬ 
out brackets) and m t (given in brackets) depending on the parameter q = yj 1 — T 2 /T c 2 . 
The dashed lines are for m t = 174 GeV, the dotted lines for m t = 150 and 200 GeV. 
The regions which mainly contribute to the baryon number violation are marked by 
solid lines. 



The erasure of the baryon asymmetry 



Fig. 5 : The ratio B 0 /Bt c as a function of the Higgs mass mu for m t = 150, 174, 
and 200 GeV. From the condition that this ratio should be at least 10' “ 5 we obtain 
an upper bound for ra# in the range 60 to 75 GeV. The same calculation without 
fermion loops leads to a qualitatively similar picture but the upper bound would be 
as low as 49 GeV. 



